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Using analytical arguments and the numerical renormalization group method we investigate the 
spin-thermopower of a quantum dot in a magnetic field. In the particle-hole symmetric situation 
the temperature difference applied across the dot drives a pure spin current without accompanying 
charge current. For temperatures and fields at or above the Kondo temperature, but of the same 
order of magnitude, the spin-Seebeck coefficient is large, of the order of fcs/|e|. Via a mapping, 
we relate the spin-Seebeck coefficient to the charge-Seebeck coefficient of a negative-jy quantum 
dot where the corresponding result was recently reported by Andergassen et al. in Phys. Rev. B 
84, 241107 (201 1). For several regimes we provide simplified analytical expressions. In the Kondo 
regime, the dependence of the spin-Seebeck coefficient on the temperature and the magnetic field is 
explained in terms of the shift of the Kondo resonance due to the field and its broadening with the 
temperature and the field. We also consider the influence of breaking the particle-hole symmetry 
and show that a pure spin current can still be realized provided a suitable electric voltage is applied 
across the dot. Then, except for large asymmetries, the behavior of the spin-Seebeck coefficient 
remains similar to that found in the particle-hole symmetric point. 



I. INTRODUCTION 

Thermoelectricity is the occurrence of electric voltages 
in the presence of temperature differences or vice-versa. 
Devices based on thermoelectric phenomena can be used 
for several applications including power generation, re- 
frigeration, and temperature measurement ji Thermoelec- 
tric phenomena are also of fundamental scientific interest 
as they reveal information about a system which is un- 
available in standard charge transport measurements 
Recent reinvigoration in this field is due to a large ther- 
moelectric response found in some strongly correlated 
materials, e.g. sodium cobaltate,^ as well as due to the 
investigation of thermopower in nanoscale junctions such 
as quantum point contacts, silicon nanowires, carbon 
nanotubes, and molecular junctions.-'- Recently, ther- 
mopower of Kondo correlated quantum dots has been 
measured^ and theoretically analyzed;^ 

Importantly, the thermoelectric effects turned to be 
useful^ also for spintronics^*'. Spintronic devices exploit- 
ing the spin degree of freedom of an electron, such as 
the prototypical Datta-Das spin field-effect transistor, ^-"^ 
promise many advantages over the conventional charge- 
based electronic devices, most notably lower power con- 
sumption and heatingji^ Recently, the spin-Hall ef- 
fect was utilized to realize the spin transistor idea 
experimentally.^^ Generating spin currents plays a funda- 
mental role in driving spintronic devices. Several meth- 
ods of generating spin currents have been put forward 
such as electrical based on the tunneling from ferromag- 
netic contacts and optical based on excitation of carri- 
ers in semiconductors with circularly polarized light.— 
Another possibility is thermoelectrical injection. In a 
recent breakthrough, the spin-Seebeck effect, where the 
spin current is driven by a temperature difference across 
the sample, has been observed in a metallic ferromagnet 
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Figure I: (a) Quantum dot attached to paramagnetic leads 
held at different temperatures, (b) Energy levels in the leads 
and in the quantum dot (at the particle-hole symmetric point, 
as a function of the magnetic field). 



and suggested as a spin current source;^ 

Some of the above mentioned methods have an unde- 
sirable property that the generated spin current is accom- 
panied by an electrical current, which leads to dissipation 
and heat. In this paper we consider a device, which uses 
the spin-Seebeck effect to generate a pure spin current. It 
consists of a quantum dot in a magnetic field attached to 
paramagnetic leads, as shown in Fig. [TJa). The leads are 
held at different temperatures. In the particle-hole sym- 
metric situation, this setup generates a pure spin current 
across the quantum dot. Away from the particle-hole 
symmetric point, a pure spin current can be generated 
provided a suitable electrical voltage is applied. By tun- 
ing the temperature and the magnetic field, large values 
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of the spin-Seebeck coefficient of the order of fca/jel can 
be reached. For tlie particle-hole symmetric case this was 
first recognized in a related problem of the charge ther- 
mopower in a negative- J7 quantum dot by Andergassen 
et al. in Ref. [8]. 

The outline of the paper is as follows. First, in Sees. HTl 
and mil we present the model and introduce the notion 
of spin-thermopower and define the spin-Seebeck coef- 
ficient. In Sec. II VI we present exact results for a de- 
vice containing a noninteracting quantum dot. Then, in 
Sec. |Vl we describe the details of the numerical renor- 
malization group (NRG) calculation and, in Sec. IVII we 
present the NRG results for a device with an interacting 
quantum dot. We also derive analytical expressions for 
the spin-Seebeck coefficient in several asymptotic regimes 
and find the boundaries between those regimes. Finally, 
in Sec. IVIII we give conclusions and examine whether 
our results can be easily observed experimentally. In Ap- 
pendix A we present an intuitive explanation why a pure 
spin current is generated in the particle-hole symmetric 
point. 

The spin-Seebeck coefficient was calculated for a 
single-molecule-magnet attached to metallic electrode&ii 
and for a quantum dot in contact with ferromagnetic 
electrodeSfi^ both within the sequential tunneling ap- 
proximation. The later system was also analyzed us- 
ing the equation of motion techniqu&i& as well as the 
Hartree-Fock approximationJi All these methods fail to 
describe correctly the effects of Kondo correlations occur- 
ring at low temperatures and magnetic fields in such sys- 
tems. Thermospin effects in quasi-one dimensional quan- 
tum wires in the presence of a magnetic field and Rashba 
spin-orbit interaction were also studied theoretically 



states in the left (a = L) and the right {a = R) lead, 
which we assume to be energy independent, as appro- 
priate for wide bands in leads with constant densities of 
states. Furthermore, we will mostly consider the system 
at the particle-hole symmetric point (in experiments, the 
deviation from particle-hole symmetry is easily controlled 
by the gate voltage) where the dot level is at Cd = /J. — ^ 
in equilibrium. In what follows we set the equilibrium 
chemical potential to zero. In such a regime the spin 
of the quantum dot is quenched at temperatures (again 
we use the energy units, fcs = 1) and magnetic fields low 
compared to the Kondo temperature^l 



K 



(2) 



To describe the system away from the particle-hole sym- 
metric point we introduce the asymmetry parameter 
<5 = 



u_ 
1 



III. SPIN-SEEBECK COEFFICIENT 

In a non-equilibrium situation the distribution of 
electrons with spin a in the lead a is described 
by the Fermi-Dirac distribution function /^o. (w) = 
1/ [e^""^"")/^" -I- l] with y.aa and being the spin de- 
pendent chemical potentials and the temperature in the 
lead a, respectively [Fig. [^b)]. The electrical current of 
electrons with spin a is 



I.^- I duj [ha (^) - fRa (OJ)] Ta (w) , (3) 



II. MODEL 

Our starting point is the standard Anderson impurity 
Hamiltoniani^ 
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E 



B 



fid 



+ ^kcccl^^Ckaa + Vkacl^Ja + h.C. (1) 
kacr kaa 

Here we assume that only the highest occupied energy 
level erf in the dot is relevant for the transport proper- 
ties. Due to the electron's spin a = ±1, in external 
magnetic field this level is split by the Zeeman energy 
g^sB. Rather than by an external magnetic field, the 
dot level could also be split by a local exchange field 
due to an attached ferromagnetic electrode)^ To keep 
the notation short, in Eq. ([TJ we express the magnetic 
field in energy units, i.e., g^B = 1- U is the Coulomb 
charging energy of the dot, while F = Fq is the total 

hybridization, where —iVa — X^fc l^feal^ / ~ ^ka + iS) 
is the hybridization function of the dot level with the 



where e is the electron charge, h is the Planck constant, 
and {uj) = TT pjl^j^p-^ Aa (w) is the transmission func- 
tion of electrons with spin cr.'^^ It is calculated from the 
impurity spectral function (cj) = — -^ImGo- (uj) where 

Ga H 



w - ed - cr 2" + iF - Eg. (w) 



is the retarded impurity Green's function and the inter- 
action self-energy Eo- (uj) accounts for the many-body ef- 
fects. Introducing the average chemical potential in each 
of the leads, 

we can parametrize the temperatures and the chemical 
potentials in the leads in terms of their average values 

1 



and 
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the temperature diflFerence 

AT^Tl- Tr, 

the voltage 

eV = ^J,L- M-R: 

and the spin voltage 

Assuming AT, eV and eVg to be small, the electrical 



current I = + 



I = 



AT 



(lit + In) — + (^ot + loi) eV + 



+ 2 (I'^t - ^ai) eVs 



(4) 



and the spin current Is = — /j. 



^•4 



AT 

(Iit-^u) — + (^ot-^o;)e^^- 



+-(Zot+2:o;)eF, 



(5) 



can both be expressed in terms of the transport integrals 

= [ dc^c^"[-/'H]r,H (6) 



where / (w) is the Fermi-Dirac function at T and fi. 

At the particle-hole symmetric point one has lot ^ 
Xqi = lo and lit — ~-^H = ^1 d^is to the symmetry of 
the spectral functions, (w) = (^^)- Thus, 
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2e / AT 
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In the absence of a voltage applied across the dot only 
the spin current will flow. In Appendix A we present an 
intuitive explanation of this result. 

The spin current is thus driven by the temperature 
difference and the spin voltage. The appropriate measure 
for the spin thermopower, i.e., the ability of a device 
to convert temperature difference to spin voltage, is the 
spin-Seebeck coefficient Sg, expressed in what follows in 
units of fes/lej, given by the ratio of the two driving 
forces required for the spin current to vanish. 
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Note that an asymmetry in the coupling to the leads, 
Tl 7^ Yr, does not inffuence the value of the spin-Seebeck 
coefficient provided F = -I- Yr stays constant. 



Through a mapping that interchanges the spin and 



the 



pair degrees of freedom, d| — )■ dj, and c 
Hamiltonian ([1]) at the particle-hole symmetric point in 
a magnetic field B is isomorphic to a negative- J7 Hamil- 
tonian away from the particle-hole symmetric point and 
in the absence of the magnetic fieldj^"— The spin-down 
spectral function transforms as A^ (w) — >■ A\^ (— = 
^t (w). Consequently, the charge thermopower in the 
transformed model. 



eV 
AT 



7=0 



TIo' 



is (up to a factor of two, which could be absorbed in a 
redefinition of the spin- voltage) the same as the spin ther- 
mopower in the original model. If transformed appropri- 
ately the results for the particle-hole symmetric case pre- 
sented in this paper agree with those reported earlier for 
the charge thermopower in a negative-?/ quantum dot.« 

Away from the particle-hole symmetric point one needs 
to apply an electrical voltage across the dot in order to 
make the electrical current vanish. In this regime, the 
spin-Seebeck coefficient 



AT 



1=0,1 =0 



1 

T 



and the required electrical voltage 



eV 



AT fx 



T 



ill 

^^Dt 
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(8) 



(9) 



can be readily derived from Eqs. (j4]) and ([5]). By means 
of the above transformation, these results can also be 
applied to the case of a negative- [/ quantum dot in a finite 
magnetic field. Notice that the applied electrical voltage 
(charge-Seebeck coefficient) maps to a spin- voltage (spin- 
Seebeck coefficient) of the negative-C/ device. 



IV. NONINTERACTING QUANTUM DOT 

The spin-Seebeck coefficient of a noninteracting quan- 
tum dot, [/ = 0, can be related to the Seebeck coefficient 
of a spinless problem with the impurity Green's function 
G (oj) = l/{uj — ed + *r). The corresponding spectral 
function is of a Lorentzian form. The Seebeck coefficient 
can be expressed in terms of the trigamma functioit^^ 



S (erf) = 27r 



Im{ 



2itT 



ItvT li 



(10) 



and is an odd function of e^, S (— ed) — —S {td)- 

Below we first study the particle-hole symmetric case 
and then generalize the results to the asymmetric prob- 
lem. 
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Figure 2: Temperature and magnetic field dependence of tlie 
spin-Seebeck coefficient (in units of fc_g/ |e|) of a noninteract- 
ing quantum dot. Tlie majcima at B = F, T = 0.34F and 
T = 0.345 are indicated witii daslied lines, (a) Exact result, 
Eq. (fTTj) . (b) Interpolation formula, Eq. (fTT)) . 



Particle-hole symmetric point, S = 

In the absence of interaction the impurity Green's func- 
tion is Gcr (lu) = l/(uj — aB + iT). We introduced the 
energy level shift B due to the magnetic field. In the 
non-interacting case, B — ^. From Eqs. ([7]) and pop 



Ss = S{i3) - S{-B) = 2S{b) 



(11) 



For T <^ max ^^F, the spin-Seebeck coeflicient is 
proportional to the temperature. 



47r2 BT 



(12) 



a results which can also be derived by performing 
the Sommerfeld expansion of the transport integrals in 
Eq. ([5]). It increases linearly with the field in low mag- 
netic fields, B <C F, 



S, = 



47r^ BT 

3 r2 ' 



(13) 



and is inversely proportional to the field in high magnetic 
fields, B > F, 



Ss 



4^r 

~B' 



(14) 



It reaches its maximal value at i? = F. 

For B ^ max (F, T) the spin-Seebeck coefficient is pro- 
portional to the magnetic field: 



(15) 



In the low temperature limit, T <C F, where we re- 
cover Eq. p3p . it increases linearly with the temperature. 



reaching its maximal value at T = 0.34F. It is inversely 
proportional to the temperature in the T 3> F limit: 



S 



(16) 



As shown in Fig. [5] an approximation that smoothly 
interpolates between the regimes of Eqs. (|T3)) . (|T4)) and 



4tt^ BT 




S7' = 



is qualitatively correct also at B,T ^ F. According to 
this interpolation formula, the spin-Seebeck coefficient at 
B = T, T = 0.34F where the maxima merge is Sg = 1.45, 
which compares well with the exact result of Eq. (|lip . 
Ss ~ 1.61. From this point, the high spin-Seebeck coef- 
ficient ridge of Ss ^ I continues to higher temperatures 
and fields along the T = 0.34i? line. Note also that at 
the ridge the magnetic field and the temperature are of 
the same order of magnitude. 



Away from the particle-hole symmetric point, ^ 7^ 

Here the spin-Seebeck coeflicient, Eq. (|8]), 

Ss^S{S + ]3) -S{S-I3), 

and the electrical voltage required for the electrical cur- 
rent to vanish, Eq. ([9]), 



eV 
AT 



S{B + 6)-S{B~d), 



can also be related to the Seebeck coefficient of the spin- 
less problem, Eq. (fTO|) . 

In Figs.[3Ja), (c) and (e) the spin-Seebeck coefficient Ss 
is plotted for various values of the asymmetry parameter 
S. For (5 < F the spin-Seebeck coefficient qualitatively 
resembles its (5 = behavior. At (5 = F it is suppressed 
in the i?, T <C F region. For (5 > F the spin-Seebeck coef- 
ficient becomes negative at low temperatures and fields. 
Its sign changes at B,T ^ y/S^ — T^, as one can check 
by requiring S[S + B) = 5(5 — B) and using the ap- 
proximate expression of Eq. ([T7]) together with Eq. PT|) . 
5 (erf) = i 5",!^^^^. 

As evident from Figs. |3jb) (d) and (f), the electrical 
voltage required to stop the electrical current, —eV/AT, 
can take both positive and negative values. At a given 
temperature this voltage is zero at a magnetic field B 
where S[B + 6) = S[B — 6). This condition is ful- 
filled close to the point where the spin-Seebeck coeflicient 
reaches a maximum as a function of the field, as shown 
with dashed lines in Figs. [S^a), (c) and (e). 
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Figure 3: Temperature and magnetic field dependence of tlie 
spin-Seebeck coefficient (in units of fc_g/ |e|) of a noninteract- 
ing quantum dot away from the particle-hole symmetric point 
for (a) 5 = r/2, (c) 5 = r and (e) 5 = 2V. The corresponding 
electrical voltages, — eV/AT, required for electrical current to 
vanish are shown in (b), (d) and (f). Along the dashed lines 
in (a), (c) and (e) the required electrical voltage is zero. 



METHOD 



To evaluate the transport integrals, Eq. ([6]), in the in- 
teracting case we employed the numerical renormaliza- 
tion group (NRG) method—"—. This method allows to 
compute the dynamical properties of quantum impurity 
models in a reliable and rather accurate way. The ap- 
proach is based on the discretization of the continuum 
of bath states, transformation to a linear tight-binding 
chain Hamiltonian, and iterative diagonalization of this 
discretized representation of the original problem. 

The Lehmann representation of the impurity spectral 



function is — 



e T 



X \{p\di\r)\^ 5 {lj ~ [Ep ~ Er)) . 

Here p and r index the many-particle levels \p) and \r) 
with energies Ep and E^., respectively. Z is the partition 
function 



This can also be expressed ai 



,30.31 



X \{p\dl\r)\''5{^-{E,-E,)), 

where /(w) is the Fermi-Dirac function. The transport 
integrals ([6]) are then calculated as^i^ 



(18) 



e T -\- e T 



In this approach it is thus not necessary to calculate the 
spectral function itself, thus the difficulties with the spec- 
tral function broadening and numerical integration do 
not arise. Furthermore, a single NRG run is sufficient 
to obtain the transport integrals in the full temperature 
interval. The calculations have been performed using the 
discretization parameter A = 2, the truncation cutoff set 
at 12a; AT, and twist averaging over Nz = 8 interleaved dis- 
cretization meshes. The spin U(l) and the isospin SU(2) 
symmetries have been used to simplify the calculations. 

An alternative approach for computing the transport 
integrals consists in calculating the spectral functions us- 
ing the density-matrix NRG method'^^ or its improve- 
ments, the complete- Fock-space NRG^^"— or the full- 
density-matrix (FDM) NRG.^^ In this case, a separate 
calculation has to be performed for each temperature 
T. The approaches based on the reduced density matrix 
are required to correctly determine the high- energy parts 
of the spectral function in the presence of the magnetic 
field. "^^ Nevertheless, since the main contribution to the 
transport integrals comes from the spectral weight on the 
low-energy scale of T, which is well approximated even 
in the traditional approach, ^^"''^ there is little difference 
in the results obtained either way. We have explicitly cal- 
culated the spin thermopower for a set of parameters B 
and T using the FDM NRG and compared them against 
the results based on Eq. (|18p : the difference is smaller 
than the linewidth in the plot (a few percent at most), 
see Figs.EIb) andE^d). The FDM NRG is significantly 
slower than our approach because it requires one calcu- 
lation for each value of T and, furthermore, because each 
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Figure 4: Temperature and magnetic field dependence of 
tlie spin-Seebeck coefficient (in units of fcs/ |e|) as calculated 
with the NRG method. The positions of the maxima at 
B = ^Tk, T = 0.26Ta', and T = 0.29B are indicated 
with dashed lines. 

such calculation is significantly slower than a single cal- 
culation based on Eq. (|18p . The trade-ofF between an 
additional error of a few percent and an improvement in 
the calculation efficiency of more than two orders of the 
magnitude is well justified. 



VI. RESULTS 

We now turn to the numerical results for an interacting 
quantum dot. We choose a strongly correlated regime 
with U = sr. The Kondo temperature, Eq. (j^l, at the 
particle-hole symmetric point is then Tk = O.llF. 

In Fig. m the spin-Seebeck coefficient is plotted for a 
range of temperatures and magnetic fields covering all 
the regimes of the symmetric Anderson model. It is evi- 
dent from the Figure that for temperatures and magnetic 
fields at or above the Kondo temperature, but of roughly 
the same size, there is a ridge of high spin-Seebeck coef- 
ficient (of the order of fc^/ |e|). At the point where the 
temperature and the field are both of the order of the 
Kondo temperature this ridge splits into two: one at the 
temperature of the order of the Kondo temperature with 
the spin-Seebeck coefficient decaying linearly as the field 
vanishes and the other at the magnetic field of the order 
of the Kondo temperature, again decaying linearly as the 
temperature vanishes. The ridges separate three slopes 
where the spin-Seebeck coefficient vanishes in directions 
away from the ridges. 




Figure 5: (a) Magnetic field dependence of the spin-Seebeck 
coefficient for several values of the temperature. The low 
magnetic field approximations are shown with dashed lines 
for T = O.OOITk, T = OMTk, T = O.ITk and T ^ Tk 
[Eqs. ini and ^\ and with long dashed lines for T = lOTx 
and T = IOOTa- [Eq. (gHU- The high magnetic field asymp- 
totics, Eq. (|25p . are shown with dotted lines, (c) Tempera- 
ture dependence for several values of the magnetic field. The 
low-temperature approximations are shown with dashed lines 
for B = O.OOITa', B = O.OITa, B = O.lTx, B = Tk and 
B = IOTa [Eqs. JTTI and and with a long dashed line 

fori? = IOOTa [Eq. (|25|) [. The high-temperature asymptotics, 
Eq. (|28p . are shown with dotted lines. In (b) and (d) the spin- 
Seebeck coefficient is presented in linear scale for T = Tk and 
B — Tk, respectively. Black dots show the results obtained 
using the FDM NRG method. 
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We show that a treatment similar to that for a nonin- 
teracting quantum dot in Sec. lIVl is possible in the Kondo 
regime (and its vicinity) provided that we account for the 
temperature and field dependence of the position and the 
width of the Kondo resonance. On top of that we find 
simple asymptotic expressions for the spin thermopower 
far away from the Kondo regime. In Fig. Elja) we com- 
pare these analytical results to the NRG data for various 
temperatures as a function of magnetic field on a log- 
log plot. In Fig. [5]Jb) one set of data is reproduced on 
a log-linear plot. In Figs, ^c) and (d) we present such 
comparisons for various values of the magnetic field as 
a function of temperature. In the regions of their valid- 
ity, the analytical expressions are found to reproduce the 
exact values reasonably well. 



Edcr = cr-ji?, where the Wilson ratio is i? = 2 due to 
residual quasiparticle interaction^ while its half- width is 
of the order of the Kondo temperature, Fo- = ^Tki^^^ 
Thus the spin-Seebeck coefficient in the B <C Tk and 
T <C Tk regime increases linearly with both the temper- 
ature and the magnetic field: 



12 r| 



(20) 



In the intermediate regime, Tk B <^ U , the mag- 
netic moment is localized but, due to strong logarith- 
mic corrections, the field dependence of and is 
nontrivial4^ Anyway, a good agreement with numerical 
data is obtained for i? < F by allowing for the broadening 
of the quasiparticle level due to the magnetic field,— 



Particle-hole symmetric point, 8 = 

We first study low temperature T <^ max (T/^ i?) 
Fermi-liquid regimes and low magnetic field B <C 
max (Ta'jT) regimes of the particle-hole symmetric 
model separately. Then we combine the results and 
present a unified theory which is capable of describing 
the spin thermopower at B,T ^ Tk- 



-T 



K 



B\ 



(21) 



and keeping the zero- field expression for the level shift, 
Ida = crB: 



BT 



3 2B'^ + (^Tk)'' 



(22) 



Low-temperature regimes, T <C max {Tk,B) 

Due to many-body effects there is an energy scale in- 
duced in the Anderson model which plays the role of 
an effective bandwidth. This energy scale is the Kondo 
temperature Tk at low magnetic fields B <^ Tk and 
rises up to F as the magnetic field kills the Kondo ef- 
fect. At temperatures low compared to this energy scale, 
i.e., in the Fermi-liquid regime, only the spectral func- 
tion in the vicinity of the Fermi level is relevant to the 
calculation of the transport properties. Here the Green's 
function can be parametrized in terms of the Fermi-liquid 
quasiparticle parameters, (w) = z^jiio — Sda + i^a), 
where Zu = [1 — (m)] ^ ^^e wavefunction renormal- 
ization factor, while Eda ~ Za [e^ -I- cr-f- + Scr (/i)] and 

Fcr = ZaT are the quasiparticle energy level and its half- 
width, respectively.^^ This parametrization provides an 
accurate description of the behavior of the system in the 
vicinity of the Fermi level, thus it is suitable for study- 
ing the transport properties; however, the features in the 
spectral function away from the Fermi level are not well 
described. 

The spin-Seebeck coefficient is derived by performing 
the Sommerfeld expansion of the transport integrals in 
Eq. ([7]) , resulting in the spin analog to the Mott formula. 



Ss 



3 At(0) 



£2 



F^ 



(19) 



In the strong coupling regime, B ^ Tk, the quasi- 
particle level is shifted away from the chemical potential. 



This expression reaches its maximal value at B = ^^Tk. 
In the Tk <C S <C F regime the spin-Seebeck coefficient 
is inversely proportional to the magnetic field: 



Ss = 



2_^T 
~~B' 



(23) 



In the high magnetic field regime, B > U, the Kondo 
resonance is no longer present and the Hubbard I approx- 
imation to the Green's function may be used,— 



(^) 



1 - (nd-a) 



(nd-a) 



-ed-U -a§+iT 
(24) 

where (nda) is the occupation of the dot level with spin a 
electrons. In this regime, {udi) and (n^^) approach 1 and 
0, respectively, and ida « cr -I- -f ) , Fo- = F. The spin- 
Seebeck coefficient increases linearly with temperature, 
but is inversely proportional to the magnetic field for 
S > J7: 



Ss 



St:' T 
3 U + B' 



(25) 



In Fig. Owe compare this expression to the NRG data and 
find an excellent agreement in the appropriate regime. 



Low magnetic field regimes, B <^ max {Tk, T) 

At low and intermediate temperatures, T ^ F, only 
the Kondo resonance will play a role in determining the 
value of the spin-Seebeck coefficient. The effective dot 
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level width, being equal to ^Tk at low temperatures. 
T <C Tfc, increases with temperature)^!^ 

2 



r 



It. 

TT 



+ W 



(26) 



Putting this expression, together with the low tempera- 
ture expression for the level shift B = B, into the non- 
interacting quantum dot formula Eq. (|15p . we recover 
the Fermi-liquid asymptotic result Eq. ([^0]) for T ^ Tk, 
while for T 3> Tk, as the width of the dot level becomes 
proportional to temperature, an asymptotic formula dif- 
ferent from that of a noninteracting level, Eq. p6p . re- 
sults: 



S, = 0.54—. 



(27) 



Between these asymptotic regimes the spin-Seebeck co- 
efficient reaches its maximal value at T = Q.26Tk- 

In the free orbital regime, T > U, the Hubbard I ap- 
proximation, Eq. ([M)) . can be used again. Now, due to 
high temperature {ridi) = {nd\) ~ i and the effect of the 
magnetic field is only to shift the spectral functions. 



A, (w)| 



B 



A„ 



B 



5=0 



As — /' (a;) sa ^ in the region where the spectral density 
is appreciable, we get a simple expression 



2 / Aujuj A-^ (i 





_ B 






-f) 



IB=0 



Lb=o 



Taking into account that in the absence of the magnetic 
field the spectral function is even, one can readily show 
that the spin-Seebeck coefficient is given by 



(28) 



thus recovering the non-interacting expression Eq. 1161 In 
Fig. [S] we show that this expression correctly describes 
the asymptotic behavior of the spin-Seebeck coefficient 
in this regime. 



Unified theory for B,T <^T regimes 

Here the physics is governed by the Kondo resonance in 
the spectral function and its remnants at temperatures 
and magnetic fields above the Kondo temperature. A 
unified description of low and intermediate temperature 
and magnetic field regimes, i?, T ^ F, is obtained using 
the noninteracting quantum dot spin-Seebeck coefficient 
formula, Eq. pT|) . and taking into account that the width 
of the Kondo resonance depends on both the temperature 
and the field. The appropriate phenomenological expres- 
sion is a generalization of the low temperature, Eq. (j26|) . 
and low magnetic field, Eq. (|2ip . widths: 



rt = ( -Tk 

TT 



^2 + (ttT)^ 



(29) 



Again we use B = B. As evident in Fig. [51 this gives 
quite an accurate approximation, reproducing the correct 
position, width and, to a lesser extent, height of the peak 
in the spin-Seebeck coefficient. It becomes even more 
accurate for U — 16T (not shown here) where the Kondo 
energy scale is better separated from higher energy scales, 
resulting in a better agreement with the NRG data in the 
B < T/f and < T < F regime. 

At the point where the B = ^Tk and T = 0.26Tk 
maxima merge, the spin-Seebeck coefficient reaches a 
value of Ss = 0.53, while the current approximation gives 
Ss = 0.63. This universal value (provided the Kondo en- 
ergy scale is well separated from higher energy scales, 
Tk <C F) is somewhat lower than that of a noninteract- 
ing quantum dot which is a direct consequence of the 
widening of the Kondo energy level. From the merging 
point the spin-Seebeck coefficient increases only slightly, 
Ss ^ 1, along the T = 0.29B line. 

As in the noninteracting case, Eq. (|T7| . we can provide 
an analytical approximation for the spin-Seebeck coeffi- 
cient by interpolating between the three asymptotic ex- 
pressions of Eqs. (HOI), (123 and (H?]): 



S7' 



TT-* BT 
12 71 



2^T 
~B 



0.54- 



(30) 



Away from the particle-hole symmetric point, S ^ 

In Fig. [6] we present the spin-Seebeck coefficient for 
three different values of the asymmetry parameter S rang- 
ing from the Kondo to the mixed valence regime. In the 
Kondo regime, S = T = ^ in Fig. (SJa) and S = 2T = ^ 
in Fig. inUb), the Kondo peak in the spectral function is 
pinned to the chemical potential in the absence of the 
magnetic field, e^^ ^ Fg..''* The Kondo temperature in- 
creases with 6, Tk oc exp (ttS^ /2UT) ^ which causes the 
positions of maxima in the spin-Seebeck coefficient to 
shift to higher temperatures and fields. Provided we take 
this effect into account, the behavior of the spin-Seebeck 
coefficient reproduces that of the particle-hole symmet- 
ric situation in Fig. 3] In the mixed valence regime, 
(5 = 4F = -J in Fig. [6{e), the Kondo peak has merged 
with the atomic peak resulting in a resonance in the spec- 
tral function of width Tg- ^ T at ida ~ Fo- at B = 04^ 
The maxima shift to T, i? ^ F and we reproduce the 
characteristic suppression of the spin-Seebeck coefficient 
at low temperatures and fields observed in the noninter- 
acting case in Fig.[3fc). In the empty orbital regime (not 
shown) the peak of width F would shift even further away 
from the chemical potential and we would reproduce the 
noninteracting result of Fig. ^e). 

The electrical voltage required to stop the electrical 
current matches the prediction of the noninteracting the- 
ory in the mixed valence regime. Fig. ^f), while in the 
Kondo regime, Figs.lHIb) and (d), it changes sign twice as 
a function of temperature at low magnetic fields. Such a 
behavior is consistent with the temperature dependence 
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I ^^"^/^ 


\ / 







BITk 



10-2 JQ-I ^qO jq2 10-1 jqO jqI jq2 

(f) 10^- 




10-2 JQ-I ^qO ^q2 jq-2 jq-1 jqO jqI jq2 

Figure 6: Temperature and magnetic field dependence of the 
spin-Seebeck coefficient (in units of ks / |e|) as calculated with 
the NRG method away from the particle-hole symmetric point 
for (a) 5 = r, (c) 5 = 2r and (e) 5 = 4r. The corresponding 
electrical voltages, —eV/IST, required for electrical current 
to vanish are shown in (b), (d) and (f). Temperatures and 
magnetic fields are shown in units of the Kondo temperature 
in the particle-hole symmetric point. Note that because of 
the increase of the Kondo temperature with 5, the maxima 
are shifted to higher temperature and magnetic field values. 
The low temperature and magnetic field structure in (e) is an 
artifact of the numerical method. 



of the charge thermopower of a quantum dot studied by 
Costi and Zlatic in Ref. Q. 



VII. SUMMARY 

We demonstrated that the spin-Seebeck effect in a sys- 
tem composed of a quantum dot in a magnetic field, 
attached to paramagnetic leads, can be utilized to pro- 
vide a pure spin current for spintronic applications pro- 
vided the dot is held at the particle-hole symmetric 
point. By tuning the temperature and the field, the spin- 



Seebeck coefficient can reach large values of the order of 
kB/\e\. Namely, for temperatures higher than Q.2QTk, 
such a large spin thermopower is available at the mag- 
netic field where T = 0.295. We carefully analyzed 
different regimes and derived analytical formulae which 
successfully reproduce and explain the temperature and 
magnetic field dependence of the spin-Seebeck coefficient 
calculated numerically with NRG. 

Replacing the magnetic field with the gate voltage, 
the same conclusions are also valid for the charge ther- 
mopower of a negative- J7 quantum dot, for which the 
charge-Seebeck coefficient of the order of ks/ \ e\ was re- 
cently reported in Ref. [8j . 

We also studied the spin-Seebeck coefficient away from 
the particle-hole symmetric point. Our analysis shows 
that in the Kondo regime the results do not change qual- 
itatively provided we take into account the increase of 
the Kondo temperature and apply a suitable electrical 
voltage across the dot to stop the electrical current. 

As the Kondo temperature in quantum dots can be 
made quite low, the magnetic fields required to reach the 
large spin thermopower regime should be easily achiev- 
able in experiment. 
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Appendix A 

According to Eq. ^ the electrical current of spin a 
electrons is determined by the energy distribution of in- 
coming electrons fL^^uj) — fji„ (lo) and the probability 
that those electrons are transmitted through the quan- 
tum dot Ta (w). In the presence of a magnetic field, 
Tcr {'^) is asymmetric about the chemical potential. As- 
suming the particle-hole symmetry and i? > it is larger 
for spin-up electrons immediately above the chemical po- 
tential (cj > 0) than for those immediately below the 
chemical potential < 0) , and vice versa for spin-down 
electrons as Ti (w) = Tf (— t^)- 

In the linear response regime we can study the effects 
of AT, eV and eVs separately as follows. 

In the presence of a temperature difference AT, 
Fig. [TJa), the energy distribution of the incoming elec- 
trons is the same for spin-up and spin-down electrons. 
For AT > the incoming electrons with a; > originate 
from the left (hot) lead while those with a; < originate 
from the right (cold) lead. Now for w > there will be a 
surplus of spin-up electrons reaching the right lead due to 
the asymmetry of 7^ (w) while for w < the same surplus 
of spin-down electrons will reach the left lead, resulting 
in a zero electrical current and a finite spin current across 
the dot. 
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Figure 7: Currents generated by (a) temperature diflerence, 
(b) electrical voltage and (c) spin voltage applied across a 
quantum dot in a magnetic field B > 0. The currents are 
shown separately for spin-up (violet arrows) and spin-down 
electrons (green arrows), as well as for electrons with energies 
higher (upper set of arrows) and lower (lower set of arrows) 
than the equilibrium chemical potential (black dotted line). 



In the presence of an electrical voltage eV > 0, 
Fig. El^b), all the incoming electrons originate from the 
left lead, their energy distribution is an even function of 
cj and is again the same for both spin projections. For 
cj > a surplus of spin-up electrons is reaching the right 
lead while for co < the same surplus of spin-down elec- 
trons is reaching the right lead. The result is a finite 
electrical current and a zero spin current. 



In the presence of a spin voltage eVg > 0, Fig. [TJc), 
spin-up electrons originate from the left lead and spin- 
down electrons originate from the right lead. Energy dis- 
tributions of the two species of incoming electrons, like 
the transmission probabilities 7^ (u), are related by re- 
flection symmetry with respect to w = 0. Consequently, 
the number of spin-up electrons reaching the right lead is 
the same as the number of spin-down electrons reaching 
the left lead, i.e. a zero electrical current and a finite 
spin current. 
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